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ORDER BOOK 
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Abstract. We introduce and treat rigorously a new multi- agent model of the 
limit order book. Our model is designed to explain a behavior of the market when 
new information affecting the market arrives. Our model has two major features. 
First, it constitutes a nonlinear Markov process. Second, it has two additional 
parameters which we call slow parameters. These parameters measure mood of 
two groups of investors, namely, bulls and bears. We explain the intuition behind 
the equations and present numerical simulations which show that behavior of the 
model is similar to the behavior of the real market. 



Starting from Louis Jean-Baptiste Alphonse Bachelier people tried to model be- 
havior of the market using various stochastic processes. Bachelier himself in 1900 
used a Brownian motion, [3J. Later people began to use Markovian diffusions, 
diffusions with jumps and even processes determined by infinite divisible distri- 
butions. Models of the price change are numerous as well as the behavior of the 
market which changes its statistics over time. 

In the last 15-20 years all US equity and futures exchanges made the order book 
electronic and information about active orders in the book is available to market 
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1. Introduction. COB and volume. 
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participants. The term market microstructure was born though the definition 
varies depending on the author, [7]. 

In the recent years numerous models of the market microstructure were intro- 
duced and studied. Here we have to mention papers [HI [19] in financial literature 
on various models of the book and limit order markets. At the same time in 
physics literature can be found a large number of the so called multi-agent models, 
see reviews [21 [121 D2] and original paper [T7]. In these models idealized market 
participants or agents submit buy or sell orders (limit or market) to the matching 
mechanism and these orders are executed according rules of exchange. Since limit 
orders are waiting for execution in FIFO queues these multi-agent models should 
be treated as a part of queuing theory. Also we have to mention the papers of 
[U EJJ where stochastic models of the order book dynamics are considered. 

In this paper we introduce and study a new multiagent nonlinear Markov model 
of the order book. Our model is designed to simulate a real life phenomena well 
known to traders in equity and futures market. The phenomena is depicted on 
the Fig 1 which contains a graph of the price of the index S&P500 future contract 
ESM08 (HLOC one minute bars) and another graph of trading volume on Friday 
April 4 of the year 2008. This is the first Friday of the month and information 
about employment is released at 8:30 am. 
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Fig 1. 

In absence of the news before 8:30 the price changes in a diffusive manner and 
volume of trading is fairly low. When news hit the market the price changes by a 
jump and then slowly returns to the initial range. At the beginning of this process 
the volume of trading increases significantly. This can be seen from the second 
graph. When time goes by the volume decreases and the price returns to the 
vicinity of the price before the announcement. Such mean reverting behavior is 
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well known to intraday traders who follow the standard calendar of announcements 
of economic indicators, [16J . 

The price change in our new nonlinear Markov multiagent model is produced 
by interaction of market participants with different roles through the matching 
mechanism. This stochastic process is a nonlinear discrete analog of the classical 
Ornstein-Uhlenbeck process. Our approach is conceptually similar to the attempt 
to explain classical Brownian motion through mechanical model of a heavy particle 
interacting with the ideal gas, see [20] . 

Our model has two major features. First, it uses ideas of kinetic theory. Specif- 
ically the ideas from the theory of Vlasov equation. Originally this equation was 
written for plasma where ions interact with long range Coulomb forces. Therefore 
an interaction between ions can not be neglected. The force acting on an ion can 
be computed by averaging the potential over the distribution of other ions in the 
configuration space. A mathematical model of this phenomena was introduced by 
H.McKean in the form of nonlinear Markov process, |15j . In our model the dynam- 
ics of parameters at any moment of time is also determined through the quantities 
which are obtained by averaging over the distribution of price at this moment. 

Second, the important feature of our model is that it introduces two "slow" 
parameters L(t) and M(t) which are functions of time. It is known that there 
are various time scales in time series or data generated by financial markets. For 
concise review see Functions L(t) and M(t) measure sentiment of the market 
(bulls and bears) about fair price of the security. They change relatively slow 
compare to the price which is the fast parameter. 

The numerical simulations of our model produce the graphs of functions L(t), M(t) 
price and volume which are depicted on Fig 2. 
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2. Description of the model. 

2.1. Matching mechanism. Here we will describe matching mechanism. We 
adopt the following conventions. A queue is represented by a half shown on 

Figure 3. Orders contain integer number of elementary unites (contracts/stocks) 
and they are executed by FIFO rule. Orders can be executed partially. 



Last First 
^ | 2 | 5 | 7 | 23 | 2 \ + 

Fig 3. 

Order book is a structured list of interacting queues as shown on Figure 4. Price 
levels are indexed by integers numbers. Each price level has two queues. One 
queue contains orders to buy not higher then this price level, the second contains 
orders to sell at price not lower then the price level. 
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Fig 4. 

The book has this form due to the following rule. If some price level X has 
k orders to buy and some price level Y has p orders to sell then if X greater or 
equal to Y them min(k,p) orders are executed immediately and removed from the 
system. This is shown on Figure 5. 
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The real book on NYSE and NASDAQ shows all orders at all levels. On CME the 
book shows only aggregated levels. 



2.2. Market participants. In the real market traders submit market and limit 
orders. The behavior of market participants determine the state of the book and 
evolution of the price. 

In our model these trading activities are separated between four groups of 
traders. 

Ask maker - Ask(p,s). 

Bid maker - Bid(p,s). 

Ask taker - Buy(s). 

Bid taker - Sell(s). 

We assume that each group consists of one trader. Ask/Bid makers fill each level 
of the book with just one contract as shown on Fig. 6. 
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Fig 6. 



If Ask taker submit buy order then she buy at the price n. If Bid taker submit 
market sell order then she sells at the price n-1. If ask/bid taker send the order 
and take contract on the ask/bid then immediately bid/ask maker fill the emptied 
level. Therefore dynamics of the book is determined by ask/bid takers. 
Ask taker (buy market orders) sends the order with exponential rate 

A n = A^ = e" c(n - L) . 
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Bid takers (sell market orders) sends the orders with exponential rate 

,, - ,,M _ p c{n-M) 
F"n — H'n ~ e 

These rates produce a continuous time Markov chain X t on Z 1 . The functions L(t) 
and M(t) vary in time according to the differential equations 

^ = -(E nXt \ L Xt - E^X S X ) + ~ n), (2-1) 



The expectation E- Kx with respect to the distribution of the process X t in the 
moment of time t. The expectation E w a with respect to the equilibrium distribution 
of the process X t with s = L = M. The quantity E^ B \ X = E n sfi s x is a constant. 
The variables A's and B's are random. They arrive at random times r's and 
#'s. We can assume that these times are exponentially distributed. We call terms 
incorporated under the sum signs an exterior forces. The triple 

(X t ,L(t),M(t)) 

forms a nonlinear Markov process. Its transition probabilities depend on the dis- 
tribution TCx t and the distribution of the exterior forces. Without exterior forces it 
can be described as a nonlinear discrete analog of the classical Ornstein-Uhlenbeck 
process. 

Questions regarding proofs of existence and convergence to equilibrium are con- 
sidered in [18]. In the following sections we explain intuition behind these equations 
and some of their properties. 

Apparently behavior of idealized market participants can be defined differently. 
In some cases the book can be analyzed only numerically, [5]. 

3. Ehrenfest model and continuous OU. 

We start with the simplest discreet analog of the classical Ornstein-Uhlenbeck 
process which also exhibits mean reverting behavior. In 1907 P. and T. Ehrenfest 
[6] introduced a model which later became known as the "Ehrenfest urn model". 
Fix an integer N and imagine two urns each containing a number of balls, in a such 
way that the total number of balls in the two urns is 2N. At each moment of time 
we pick one ball at random (each with probability 1 / 2 N) and move it to the other 
urn. If Y t denote number of balls in the first urn minus N then Y t , t = 0,1, 2...; 
form a Markov chain with the state space {— N, X}. 

This Markov chain is reversible with binomial distribution as a stationary mea- 
sure 

2AH (\ 
(N + k)\(N — k)\ 

and the generator 



k = -X V: 



K E f{k) = i (l - A) f(k + 1) + 1 (l + 4 ) /(* - i) - Mh 



2 V X Jy '2V N 



In 1930 Ornstein and Uhlenbeck [10J introduced a model of Brownian particle 
moving under linear force. It is a Markov process with the state space R 1 and the 
generator 

d 2 d 



dx 2 dx ' 
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The OU process is reversible with invariant density 

tt ou (x) = ^e~ cx \ 

The two Markov processes are related. Under scaling 

k ~ x, 

spatial 1 

in 



N "'" 



i 
c 



the generator 21 s takes the form 

2 \ \/n J Jn 2 V \/n I Jn 



Using Teylor expansion after simple algebra we have 

K E f{x) = ^ [f"(x) - 2cxf\x)\ + ... = i-2t oc/ + ... 

That shows on a formal level that St^converges to the generator 2l OC7 . Also bino- 
mial distribution n§ converges to a Gaussian measure with the density 

4. Discrete OU process. Speed and jump measure. 

Now we introduce a different discretization of the classical OU process. In its 
simplest form X t is a continuous time Markov process on Z 1 . The intensity of a 
jump up is 

A = X L = e ~ c(n ~ L) 

and of a jump down is 

,. _ ,.M _ p c(n-M) 

where c > and L and M are real numbers. Later we will make L and M time 
dependent. 

Now we want to study behavior of the process specifically the probability of 
getting to infinity and expectation of time of returning from infinity. We assume 
that L = M = 0. Since both spatial infinities are identical we consider the 
process on right semi-axis. A method of study such processes by means of an 
auxiliary space(A, E) was introduced in the paper [9]. The set E consists of points 
{ x n}n=o c M >o, where 

x = — = 1, 
fJ>o 

1 

X!= x + — = 2,... 
^0 



_ fhJJq ■ ■ ■ _ m(n+l) 

AoAi . . . A n 

and 

Xqo = lim x n = oo. 

n— >oo 

This fact implies that the process is recurrent, see [1: 16. b]. 
The measure [i is defined by the rule 

fi n = n({x n }) = A ° Al ' ' ' Ara-1 = e _c ™ 2 . 

The ideal point entrance point, i.e. the expectation of time to come to 

the finite part of the phase space starting at infinity is finite. This follows from 
the estimate 

oo oo n oo 

n=0 n=0 fc=0 n=0 

Whence, X 4 is very close to a Markov chain with finite number of states. 
The detailed balance equations 

ir(n)X n = n(n + l)fi n +i, 7r(n)/4i = n{n - l)A n _i, 
are satisfied for the distribution 

vr s (n) = ^e" c (— ) 2 , s = ^±^. 

The normalization factor H = H(s, c) is given by 

- ( 5 c ) = e 9 — , - 

where Jacobi theta function, [5], is 

e(u,r) = J2e 2wivn+wiTn2 . 

It is interesting that invariant measure 7r which should depend on both parameters 
L and M depends on s only. The chain is ergodic and reversible. 
For L = M = the process X t has generator 

2l/(n) = e- cr 7(n + 1) + e cn /(n - 1) - (e" cn + e cn ) /(n). 

Let us show how on formal level 2t oc/ can be obtained from the generator of discrete 
process. If one scales 

n ~ x, 
spatial 1 

in 



c 

c ~ 

in 



then 



+ e 



X 



f(x-A=)-f(x) 



Using Taylor expansion we have 

*/(*) = - [f"(x) - 2cxf'(x)} + ... = -% ou f{x) + .... 
n n 

Apart from the factor ^ this is a generator of the classical model. Also under this 
scaling discrete Gaussian distribution 7r s (n) converges to the Gaussian distribution 
with the density tt ou (x). It is naturally to call X t a discrete Ornstein-Uhlenbeck 
process. 

In addition to s it is natural to introduce another variable d = ^ = ^~. The case 
when d = we call agreement and d ^ we call disagreement. The invariant 
measure does not depend on parameter d but the intensity of jumps does. Simple 
arguments show that at the equilibrium for any d 

E w sX x = E^Hx ■ 

In the case of disagreement at the equilibrium expectations are changed by expo- 
nential factor e cd 

E w sX x = e cd 'E^sX'^, E^s^i^x = G cd E^s ^i s x . (4-1) 

Apparently when d > the expectation increases and if d < decreases. 

5. Slow functions L and M. Nonlinear Markov process. 

Let us assume that parameters L and M are not constants but are some functions 
of time. They change according to differential equations 

dL(t) 



dt 



-(E nXt X L Xt - E n s\ s x ) + J2 A kS(t - Tk ), (5.1) 



The variables A's and -B's are random and positive. They arrive at random times 
r's and 0's. We can assume that these times are exponentially distributed. We 
call terms incorporated under the sum signs an exterior forces. The triple 

(X t ,L(t),M(t)) 

forms a nonlinear Markov process. Its transition probabilities depend on the dis- 
tribution Tix t °f the first component and distribution of the exterior forces. 

Lets us explain the intuition behind these equations. In reality there are two 
groups of market participants. They are called bulls and bears. Consider the bulls 
which try to make the price higher all the time and submit buy market orders 
with intensity X^. Buying a contract is an expense for them. Bulls have two 
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quantities in mind. One is the quantity E n s\ s x which is how much they are willing 
to pay per unit of time. Another quantity is the function L(t) which measures 
their expectations about fair price. If their expectations about fair price are too 
high then the quantity E nx \ Xt — E n s\ s x is positive and the bulls decrease their 
expectations according to the differential equation 

= ~( E nx t ^X t ~ E n s\ s x ). 

When positive news arrive the bulls change their expectations by a jump. That 
is contribution of the sum corresponding exterior sources. Exactly the same argu- 
ments apply to the bears which submit sell orders with intensity fi^f . 

The situation resembles random walk in random environment. The functions 
L(t) and M(t) determine transition probabilities P L - M of X t and expectations 
E L ' M with respect to these probabilities. These functions are hidden parameters 
of the system and it is possible to condition on them. Instead we condition on 
the exterior forces A's and B's which through differential equations determine 
the functions L and M. These transition probabilities are denoted by P A > B and 
corresponding expectations by E A,B . We refer to them as quenched probabilities 
and expectations. 

The Poisson measure on exterior forces A's and -B's we denote by P. The product 
measure P A ' B x P determines averaged transition probabilities of the process X t . 
The corresponding expectation we denoted by E Xt - 

6. Hydrodynamic limit. 

In this case one assumes that the variables L and M change change with macro- 
scopic time t while the process X T moves with microscopic time r = t/e where 
e > 0. In the limit e — > the process X T is always at equilibrium i.e. it has 



stationary distribution n s for all moments of time. Therefore, using 4.1 



ll 



*M = -(e cd - l)E na \ s x + J2^S(t - r k ), 

= +(e cd ~ 1)^/4 + E B ^ ~ **)■ 

Subtraction one equation from another and discarding the "news" we obtain closed 
equation for the function d(t): 

d*(t) = -{e cd - l)V, 

where V = E n s\ s x = E n sfi s x . This can be solved explicitly using substitution 
d=\\og(± V ). 

If a solution is negative for some moment of time then it stays negative for all 
moments and it is given by 

1 Ae cVt 

d{t) =c ,og J^TT- - 4>0 - 

The solution increases monotonously from negative infinity to zero when time runs 
the whole axis. 

On the opposite, if a solution is positive for some moment of time then it stays 
positive for all moments of time when its defined. The solution is given by the 
formula 

1 Ae cVt 

The solution is defined and decreases monotonously from positive infinity to zero 
for t > t = -^logA 

7. Continuum limit. 

Now we consider a continuum limit of the discrete model. The continuous time 
process X(t) with values in R 1 is defined as a solution of the stochastic differential 
equation 



dX(t) = e< h ~W dw(t) - 2c (x(t) - 



dt 



where c > and the functions L(t) and A4(t) are 

dt " 1 h 



dM{t) 



+ (e c(L-A4) _ x y 



dt 

There is a simple formal scaling relation between two models. The generator of 
the discrete model has the form 

2t/(n) = e-< n - L ^f{n + 1) + e c( "- MW) f(n - 1) - ( e - c{n - L{t)) + e c(n " MW) ) f(n) 
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Introducing new variables 



L + M , L-M 

, d= r , 



2 ' 2 
we have 

St/(n) = e cd [e- c(n - s) /(^ + 1) + e c(n ~ s) /(™ - 1) - (e~ c(n - s) + e c(n ~ s) ) /(n)] . 
If one scales 



spatial 1 



n ~ x 

1 



c 



L + M ,-fr \ a\ 

M ~ L + A/f - v^(L-.M) 



then 

s~ — - — , d~V ra (L-Aa); 
and for the generator we have 

*/(*) ~ -U< L ~ M ) [/"(*) - 2c(x - + .... 

Apart from the factor - this is the generator of the continuous model. Time 
scaling t ~ nt removes this factor. The differential equations for L and M when 
they scaled the same way produce differential equations for functions L and M. 

If L(i) = M(t) = s, then X(t) is just a classical Ornstein-Uhlenbeck process 
with mean s. It can be proved similar to [18] that if the process starts with some 
L(0) 7^ M(0), then X(t) converges to OU and L(t) and M{t) converge to the same 
constant s. 

8. Propagation of chaos and Monte-Carlo simulation. 

For the purpose numerical simulation we consider multi-particle approximation 
of the process (X t , L(t), M(t)). 

Everywhere below we will denote by c(«) any nonnegative nondecreasing func- 
tion of some parameters. Define a banach space B a consisted of all vector-functions 
P = {Pn} ne z w ith the Banach norm 



\P\ 



a 



2J \p n \ exp ( — + a\n\ J , ae 
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Also let p(0) denote the initial distribution of X t . In [T8] we have proved the 
following theorem: 

Theorem 8.1. For any initial probability measure p(0) G B a , the process X t 
is defined correctly on the half-line [0, oo) and for any t > its distribution pit) 
belongs to B a . 

To simplify notations we assume c = 1. Fix some integer N and consider 
the distribution of iV particles on : £i, £2, • • • ,£/v; £ Z. The positions 
of particles change with time & = We define a stochastic process Xn = 

(pN,Lff(t),MN(t)), where pjv(-) = PN(*,t) is empirical distribution of at the 
moment t. All particles jump with the same rates 

Av(6) = e-^~ L -\ M6)=e (6_%) - 
The functions Ljv(t), Mjy(t) change according to differential equations: 

^f= -(E„ e zPiv(n)A w (n)-i?^)+EMi-r fc ); 

The random variables A,B,t, 8's have the same distributions as before. 

We also define an auxiliary inner product on the space of infinite sequences £(n) 
and 77 (n), n G Z; as 

(f > = £( n M n ) ex P (y + a \ n \) > «eR. 

This inner product produces the euclidean norm denoted by | • | a . Now we are ready 
to prove the following theorem. 

Theorem 8.2. Let p(0) G B a+1 , where a is an arbitrary real number. Assume 
that E\\p N (0) - p(0)\\ a+1 -+0, as N ^ 00. Then for all t > 

sup £ , |pat(s) — — )■ 0, as iV — >■ 00. 

Proof. We will prove the theorem in assumption of absence of the news. For our 
purposes it will be convinient to decompose the inifitisemal matrices of any par- 
ticipants corresponding to the processes X and X^ into the sum of their diagonal 
and off-diagonal parts: 

H — H + V; H N = H N0 + V N . 

Lemma 8.3. (auxilary results) Let t > 0, then for all N G N 

(1) e L ^ < c(t); e~ M ® < c (t); e L ^ < c(f); e~ M ^ < c(t). 

(2) \e ±L - c ±Ln \ < c(t)r(t) and |e ±M - e ±Mjv | < c(t)T(t), where 

T(t)= f Y / \p(n,s)-p N (n,s)\e^ds. 
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(3) ET 2 (t)<c(t)su Ps<t E\p N -p\ 2 a . 

(4) lfp(0) e B a , then E\p(H m - H )\ 2 a < c(t) sup s < t E\p N - p\ 2 a . 

(5) Ifp(0)eB a , thenE\p(V N -V)\l<c(t)su Vs < t E\p N -p\l. 

(6) lfp(0) e B a , then E(p N - p, (p N - p)H N ) a < c(t) sup s < t E\p N - p\ 2 a . 

(7) Let E\\p N (0)\\ a+1 < oo, then 

J2^nEQ n <c(t,E\\p N (0)\\ a+1 ), 

n 

where 

Qn = PN(n - l)Aiv(n - 1) + p N {n + + 1) + p N (n)\ N (n) + p N (n)fx N (n). 



Proof. The first statement follows directly from (5.1 ) and (8.1 ). Let £ = e — e 
then (5.1) and (8.1) imply 

£ + i^A^ = J>(n) - Piv(^))e- n . 
This ODE can be solved explicitly and we get 



\e~ L - e~ Ljv | 



e~ CA * / e c ^ s ^{p{n) - p N {n))e- n ds\ < c(t)T(t), 



and therefore using the first statement 

\e L - e LN \ = e L+LN \e~ L - e~ Ljv | < c(t)T(t). 

The other statements of the second point can be obtained by analogy. Applying 
Cauchy- Schwarz inequality twice we get the third statement: 

ET 2 (t) < c(t)J*E (y2(PN-p)(n)e^A ds 

< c{t) J* E\p N - p\ 2 a e-^ + ^A ds < c(t) sup E\p N - p\ 2 a . 

We get the fourth statement in the following way: 

E\p(H m -H )\l = Ej2e^ +aH p(n) 2 \H m (n,n) - H (n,n)\ 2 < 



„2 



< Ej2^ +{a+2)H p(n) 2 (\e L -e LN \ + \e~ M - e~ MN \) 2 < 

< \J2 e ^ +ia+2)H P(n) 2 ) -c(t)ET 2 (t) < c(t) sup E\p N -p\ 2 a . 



v ngZ 

V- 



s<t 



=c(t), since ||p||a<c(t) 
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We get the fifth statement from the next estimation: 

E\p(V N - V)\l = Ej2 e ^ +aH \ e ~ n+1 ( eL - e LN )p(n - 1) + e n+1 (e~ M - e - MN )p(n + 1)| 2 < 



< const ^ +aH [e~ 2n+2 p{n - l) 2 E{e L - e L ») 2 + e 2n+2 p{n + l) 2 E{e~ M - e"^) 2 ] < 

neZ 

< const ^(e I ^ +a|n+1 ^ 2n + e^ +a|n - 1|+2n )p(n) 2 -c(t) J Er 2 (t) < c(t) sup E\p - p N \ 2 a . 

nez s ^ 



2 

a ' 



=c(t), since ||p|| Q <c(i) 

Denoting for brevity x = pn — P, we get the sixth statement 
E(x,xH N ) a = J2^ +aH Ex n (x n ^e- n+1+L * + x n+1 e n+l ~ M » - x n e~ n+LN - x n e n ~ M ^ 

n 

= J2e^ +aH - n Ee LN (ex n ^x n - x 2 n ) + ^ e f+*\n\+n Ee -M N ( eXnXn+1 _ x y < 

n n 

< c(t) e^- n + aH Ex 2 n _ x + c(t) e^- n+aH Ex 2 n+1 < c(t)E\x\ 2 a < c(t) sup E\p - p N \ 

s<t 

n n — 

Note that ^2 n a n EQ n < c(t)E\\pN(t)\\ a+ i. Fix t > and consider the Markov 
chain X on Z with intensities of jumps 

n ->• n + 1 :A n = exp {-\n\ + max (+L(0), -M(0)) + ^A^t} ■ I {n > 0} 
n n - 1 :/i n = exp{-|n| + max(+L(0), -M(0)) + E^A^t} • I {n < 0} . 

By analogy with lemma 2.4 in [18] it is easy to prove X defines a semigroup 
of transitions probabilities Px{-,-), which are bounded operators in B a+ i. Also 
due to the choice of transition probabilities we can always construct a coupling 
between one participant of Xn '■ £i and X such that the absolute value of position 
of X is not greater than \^\ for all s <t. Therefore we conclude £7||pjv(0IU+i — 
E\\p x (t)\\ a+1 = E\\p N (0)P x (0,t)\\ a+1 < c(t)E\\p N (0)\\ a+1 . □ 

Now we are ready to prove the theorem. Let us consider the process Xjy. During 
time dt the variable pjv(t) can be changed in the following way: 

{Pn(ji) - ±, with probability P out dt = Np N (n)(X N (n) + fi N (n))dt; 
Pn(ji) + with p - ty P in dt = Ndt(p N (n - l)\ N (n - 1) + Pn(ji + 1)hn(ti + 1)); 

p N (n), with p - ty 1 - P out dt - P in dt. 

The process X is deterministic, so during time dt 

p(n) ->■ p{n) + {pH{t)){n)dt = 

= p{n) + (A„_ip(n - 1) - (A n p(n) + ^ n p{n)) + ^ n +iP{n + l))dt. 
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Using that we get 



(pj\r(n) - p(n) - + 0(dt), P out c2i; 



d £| Pjv - p\l = | ^ rf +a|n| ^(Piv(n) - p(n)) 2 < 



iV 

(pjv(n) - p(n)) 2 ->• <( (p^n) - p(n) + + 0(dt), P in cft; 

(piv(n) - p(n) - {pH{t)){n)dtf, 1 - P out cft - i^dt. 
Applying Markov property and opening brackets we get 

j t E{p N {n) -p{n)f = 2E{p N {n) - p{n))(p N H N {n) - pH{n)) + ^EQ n = 

2E(p N {n) -p(n)) ■ [(p N - p)H N (n) +p(H N0 - H )(n) + p(V N - V)(n)} + ^EQ n , 

Summing up, we get 
d 2 d 

dt ElPN - pL= dt 

2E{p N - p, (p N - p)H N ) a + 

2 (E\ PN - p\l) 1/2 ■ [(E\p(V N - V)\l) 1/2 + (E\p(H m - H )\l) 1/2 ]+±J2 a " E Q- 
Summarizing results of lemma 8.3| we get the final estimate: 

j t E\p N - p\l < c(t) sup E\p N (s) - p(s)\l + l C (t). 

The Grownall's lemma implies the result. □ 

In order to produce simulations depicted on Fig. 2 we used multiparticle ap- 
proximation with iV = 1000 particles {£,i}^ =1 described above. We suppose that 
the tick is e — 0.1, so the process lives on the lattice eZ. All particles jump 
with the same rates 

where K is a regularizing parameter. We put c = 0.05 and K = 0.0333. There is 
a jump on the picture at the moment of event L — > L — A, where A = 80. After 
that at some exponentially distributed moment of time there is a second jump 
L^L+j. 
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